Generalization of the Bloch-Messiah-Zumino theorem 
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It is shown how to construct a basis in which two arbitrary complex antisymmetric matrices C and C' 
acquire simultaneously canonical forms. The present construction is not restricted by any conditions 
on properties of the C + C' matrix. Canonical bases pertaining to the generator-coordinate- method 
treatment of many-fermion systems are discussed. 
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The Bardeen-Cooper-Schrieffer (BCS) pairing theory 
0, and its generalization by Bogolyubov 0, are based 
on fermion wave functions that have the form of fermion- 
pair condensates, i.e., 



| C) = exp 
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where a+ are the fermion creation operators, |0) is the 
fermion vacuum, and C mn is an antisymmetric complex 
matrix. Up to a unitary transformation of the single- 
particle basis, 



mm' rn' ' 



state ([!]) is equal to the so-called BCS state 
\ C ) = II (l + «mc m a+a+) |0), 



m>0 



better known in its normalized form, 
\C) 
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The Bloch-Messiah-Zumino theorem provides the 
link between the two forms of state \C), Eqs. (|l|) and (j^), 
by stating that every complex antisymmetric matrix can 
be brought by a unitary transformation into its canonical 
form, i.e., 



(U CU^mn — S n C n Smni 



(6) 



where index m denotes the so-called canonical partner of 



the state m, the phase factors s*=s r 



-st have for the 



canonical partners opposite signs, and numbers c n —Cf L 
are real and positive. Standard notation m>0, used in 
Eqs. (H) and ([|), means that only one state is taken 
from each canonical pair. The proof of the theorem goes 
through a diagonalization of the hermitian matrix C + C ', 



that yields the unitary transformation U and real, non- 
negative, pairwise degenerate eigenvalues (? n . 

When using states (|l|) in applications beyond the 
mean-field approximation, and in particular in the gen- 
erator coordinate method (GCM) |^-^|, the matrix ele- 
ments and overlaps depend on the product matrix C + C . 
For example, the overlap of two states (|l|) reads ||[|||] 

(C'\C) =det 1/2 (l + C+C'), (7) 
and the transition density matrix is given by |q] 



(C'\a+a m \C) 
(C'\C) 



■1 + C+C) C+C 



(8) 



It has been realized long time ago Jl(| that the matrix 
C + C is also pairwise degenerate, which facilitates cal- 
culation of the phase of the overlap, otherwise ambigu- 
ous because of the square root appearing in Eq. (0). 
Moreover, under certain conditions it has been proved 
in Ref. O] that it is enough to give up the unitarity of 
matrix U to bring both matrices C and C simultaneously 
into the canonical forms analogous to (^). The same fact 
has later been rediscovered in Ref. [Q , although the nec- 
essary restrictions on matrices C and C have not been 
recognized there. 

In the present paper, I generalize results of Ref. JllJ by 
deriving canonical forms of two arbitrary complex matri- 
ces C and C in a common canonical basis. These results 
are not restricted by any conditions on matrices C and 

a. 

Let us begin by recalling the notion of the Jordan form 
(see e.g. Jl3|) of an arbitrary complex matrix. Focussing 
our attention on the matrix C + C", the vectors defining 
its Jordan basis can be arranged in columns of matrix 
W, and one has 
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where matrix D is block diagonal (composed of the Jor- 
dan blocks). One can attribute the number of the block 
Ii, the length of the block Li, and the number within the 
block ki, to every index i that numbers the Jordan basis 
vectors. In this notation, D has the form: 



D jt = S, , I) 



kj ki ' 



(10) 
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where in every block matrix Dj, k , reads 



(11) 



i.e., it has a common complex number Di on the main 
diagonal and the ones just above the main diagonal. 

Basis vectors belonging to a given block I form the so- 
called Jordan series of length L. The series starts with 
the basis vector called the series head, and ends with 
an eigenvector of C + C . The whole series is uniquely 
determined by the series head, because the remaining 
basis vectors in the series can be obtained by a repeated 
action of C + C on the series head. The basis vectors in a 
given series are not unique, because a linear combination 
of these vectors may give another valid series head, and 
leads to the same Jordan canonical form. Explicitly, this 
transformation reads 



W m kOLkk' 
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where the transformation matrix akk' depends on L ar- 
bitrary complex numbers otk (only a,\ must not vanish), 
and has the following explicit structure: 



Ot-kk' 
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It is easy to check that matrices having this structure 
form a group. It is also easy to check that vectors W' k 
form the Jordan series, similarly as vectors W m k do, and 
that they can replace vectors W m k in the Jordan basis, 
giving the same matrix D in Eq. (^|). 

According to the Jordan construction, the whole space 
in which acts matrix C + C splits into subspaces spanned 
by the Jordan series. The number of eigenvectors of 
C + C equals to the number of different series, or to the 
number of Jordan blocks, and is in general smaller than 
the dimension of the matrix C + C. Some matrices (her- 
mitian or not) can be fully diagonalized, i.e., they have 
numbers of eigenvectors equal to their dimensions. This 
corresponds to the case when all the Jordan series have 
the length equal 1. 

One calls two blocks degenerate, or two series degen- 
erate, if they have the same diagonal number Di, and 
they have the same length L. The latter condition is 
very important, because only degenerate series defined 
in such a way can be mixed; this is an analogue of the 
possibility to mix degenerate eigenvectors of a matrix 
which can be fully diagonalized. If two series have dif- 
ferent lengths then vectors of a longer series cannot be 
admixed to those of the shorter series, even if the series 
have the same diagonal number D/. If the matrix can 
be fully diagonalized, then all the series have length 1, 



the number of eigenvectors equals to the dimension of 
C + C", and transformation ( |l2|) reduces to the possibility 
of arbitrarily normalizing every eigenvector. 

After these necessary preliminaries, let us proceed with 
presenting the main results of the present paper. Multi- 
plying from the left-hand and right-hand sides the eigen- 
equation (|^) by W~ x , and then transposing, we obtain 
that 



C ( C+W~ 



= W- lT D T 



(14) 



which multiplied by C + from the left-hand side gives 

(C+C) (c+W~ lT ) = (c+W- lT ) D T . (15) 

One can see now that matrix C + C has another equiv- 
alent set of Jordan series, i.e., 

E( c+c ') mn v ^=J2 v ^ D ^ ( 16 ) 

n j 

where the new matrix of basis vectors V is given by 

V = C + W~ lT J. (17) 

In every Jordan block, matrix J has the ones on the 
skew diagonal and zeros otherwise, i.e., Jkk' = Sk,L-k'+i- 
Hence, when an arbitrary matrix is multiplied by J 
from the right-hand (left-hand) side, the order of its 
columns (rows) is flipped. In particular, one obtains that 
D T =JDJ. 

We can now analyze cases of different degeneracies of 
the Jordan blocks. The arguments given below closely 
follow proofs presented in Ref . Jll]] , only with the degen- 
eracies of eigenvalues replaced by the degeneracies of the 
Jordan blocks. 

Let us first suppose that C + C has a non-degenerate 
Jordan block. Then, in this block the basis vectors V 
must be connected with the basis vectors W by transfor- 
mation f[l2|) , i.e., 



( c+, ^ ,Tj :u=i> 



(18) 



fc=i 



Multiplying Eq. (|lS|) by W^ 1 from the left-hand side, 
and by J from the right-hand side, one obtains 



kk' 



i aJ )kk' 



(19) 



However, the matrix on the left-hand side of this equa- 
tion is antisymmetric, while that on the right-hand side is 
symmetric; therefore, matrix akk' must vanish. This re- 
quires that, in a non-degenerate Jordan block, all vectors 
V vanish, which contradicts Eq. ( p^ ) , unless the block has 
length L=l and Dj—0. Therefore, matrix C + C cannot 
have non-degenerate Jordan blocks apart from the sub- 
space of L=l eigenvectors with all eigenvalues equal zero. 
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One can set this subspace aside, and assume from now 
on that C + C is non-singular and has an even dimen- 
sion. In this case, matrix C + C cannot have any non- 
degenerate Jordan block, and hence Jordan blocks must 
appear in degenerate pairs. (In odd dimensions, C + C 
must have at least one null eigenvalue, which can be sep- 
arated, and the remaining matrix can be treated in the 
even dimension). 

In the present considerations, it is enough to consider 
only pairs of degenerate blocks; had the higher degen- 
eracies of the Jordan blocks occurred, one could have 
considered one pair after another, and at each step one 
could reduce the dimension of the problem. This is pos- 
sible here, and has not been possible when considering 
degenerate eigenvalues in Ref. jy], because the whole 
space can be separated into the subspaces corresponding 
to the Jordan blocks, while it cannot be separated into 
subpaces corresponding to the eigenvalues. 

Let us now consider a pair of degenerate Jordan blocks, 
each block having length L and the common diagonal ele- 
ment Dj=Dj. One can adopt here the standard notation 
that originally pertains to the canonical pairs, namely, we 
denote the indices of the two degenerate blocks by I and 
/. Similarly, indices inside these two blocks are denoted 
by k=l, 2,. . . ,L and k=l, 2,. . . ,L, respectively. Note 
that vectors in these two blocks form series, i.e., they are 
arranged in a specific order; therefore a vector at a given 
position must be associated with the vector at the same 
position in the second block. 

Since for matrix C + C two equivalent Jordan bases 
exist, W and V, vectors in series V must be linear com- 
binations of those in series W. In the pair of degenerate 
Jordan blocks, this leads to the following relations be- 
tween the two series: 

L L 

(c+W- lT j) =Y^W mk a kk , +Y j W ml l3 kk , , (20a) 

V / mk' 

k=1 k—1 
L L 

( C+W ~ lTj ) mk , = J2 W ™ k ^kk' +J2 W mk e kk', (20b) 

mk ' fc=l fc=l 

All the four matrices a, (3, 7, and e have the same 
structure (|13|) . One may now proceed with multiply- 
ing Eqs. ( |20a| ) and ( |20b| ) from the left-hand side either 
by or by W£~ , and from the right-hand side by J. 
Since all matrices aJ, (3 J, "/J, and e J are symmetric, one 
then obtains that a=e=0 and 7=—/?. 

Therefore, the canonical form of the C + matrix reads 

(w^c+w- 17 ) = S /,c£^w ( 21 ) 

where the symmetric matrix C I —C I occupies the off- 
diagonal part in every pair of the degenerate Jordan 
blocks, and 

C 1 = sj/3* J, (22) 



for j3 having form ([l3|). Following the standard nota- 
tion, we have defined the phase factors s/=— sj in such a 
way that the skew-diagonal matrix elements of C 1 (that 
are all equal one to another) are real and positive, i.e., 

qL-*+i>°- 

Since the canonical basis of C + is the same as the 
Jordan basis of C + C ', matrix C' must in the very same 
basis assume an analogous canonical form: 

(W T C'W). i = s* It C' k %5 IiIi , (23) 

where the symmetric matrix C a =C reads 

C' 1 = s*J[3'* (24) 

and /3' has also form (|r|). Finally, in order to satisfy 
Eq. (||) matrices (3 and /?' must obey the following con- 
dition: 

PP* = D. (25) 

This leaves us still some freedom in the choice of the 
canonical basis, because any solution of Eq. (^5|) gives 
one valid canonical form. Two obvious choices are, for 
example, /3=D and /?'* =1 or l3=VD and /3'* = v^D, where 
any one of the possible branches of the matrix square root 
can be taken. 

Equations ( ^l|) and (|2^) complete the proof of the 
canonical forms of two arbitrary complex antisymmet- 
ric matrices C and C . Both these matrices can be si- 
multaneously transformed by matrix W (in general non- 
unitary) into the block-diagonal forms with non-zero el- 
ements only between pairs of degenerate Jordan blocks. 

Needless to say, whenever matrix C + C can be fully 
diagonalized, which was the case in Ref. both matri- 
ces C and C acquire in the canonical basis the standard 
canonical forms analogous to Eq. (||), i.e., 

(w^C+W- 1 ^^ = Sj c*5 ri , (26) 

and 

(W T C'W) .. = sf^S fi , (27) 

where c*c^=Dj. 

In Ref. [O] it was noticed that an incorrect conjecture 
was formulated in Ref. [ fLl) , namely, the conjecture that 
the simple forms of Eqs. d26| ) and (|2^) can always be 
achieved. In the present study we have seen that these 
simple forms occur only when matrix C + C can be fully 
diagonalized. In fact, this is the case which occurs most 
often in applications. Therefore, let us now discuss con- 
ditions for the full diagonalization of C + C . 

In the applications given in Ref. , the full diagonal- 
ization of matrix C + C was secured by using a model in 
which matrices C were time-even, 

C+ = U T C T U%, (28) 
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and the Hermitian and time-even matrices C defined by 

C = -U T C (29) 

were positive definite. In these equations, Ut is a uni- 
tary and antisymmetric matrix, U^ =U T ~ 1 =—U T : - The 
positive definiteness of matrices C was in [ |i"t| guaran- 
teed by a special form of matrices C. In that study, the 
GCM states were constructed within the SCEM model 
|^5| , fl6l| , and therefore matrices C had the form shown in 
Eq. (2.13) of Therefore, the corresponding C ma- 

trices were all equal to exponents of hermitian matrices, 
and hence trivially positive definite. 

In the general presentation of the present paper, con- 
ditions ( |28| ) and ( p9| ) can be formulated as follows: If 
there exist a unitary antisymmetric matrix Ut such that 
Eq. ( pi| ) holds for C and C, and at the same time at 
least C or C is positive-definite, then matrix C + C can 
be fully diagonalized, and the simple canonical forms j26| ) 
and ( p7| ) exist. The proof of this statement has been given 
in Rcf. jll] (Appendix C), and will not be repeated here. 

The positive definiteness of C is a required condition, 
because the hermitian square-root of C must exist. Un- 
fortunately, this condition cannot be released, i.e., if both 
C and C are not positive definite, it may happen that 
matrix C + C cannot be fully diagonalized. An example 
of such a situation is provided by the following two 4x4 
matrices: 

C =(-A? 0) and C '=(-A-f)> ( 3 °) 
where the two-dimensional matrices A and A' read 

A=(y Q ) and A'=(°J). (31) 

For the standard time-reversal matrix Ut given by 

Ut = ( _° 7 J ) one has C = ( A * ° \ , (32) 

and C' has the same form. Neither C (for a^O) nor C' 
is positive definite, and the C + C'=CC matrix, 

cannot be fully diagonalized, unless a^a* . 

However, for any small but non-zero imaginary part 
of a, matrix C + C can be fully diagonalized. Therefore, 
this example also shows that the positive definiteness of 
(time-even) matrices C or C is only a sufficient condi- 
tion for the full diagonalization of C + C'=CC' , but it is 
not necessary. Moreover, it is clear that matrix C + C 
cannot be diagonalized for a=o*, because in the limit of 
Qa — >0 two eigenvectors of C + C become parallel. This 
illustrates the difficulty of diagonalizing C + C numeri- 
cally for small values of Qa; the task is then bound to 
become ill-conditioned. 



In the GCM, matrices C are most often obtained 
from solutions of the Hartree-Fock-Bogoliubov (HFB) or 
Hartree-Fock+BCS [S equations for time-even states. In 
these cases, matrices C are diagonal in the HFB or BCS 
canonical bases j§] (composed of pairs of time-reversed 
states), and their eigenvalues are equal to v m /u m =c m , 
where v m and the standard quasiparticle ampli- 

tudes of Eq. (g). Here, the canonical pairs are defined by 
the time reversal, and therefore the eigenvalues c m , can 
have, in principle, arbitrary signs. 

However, in the BCS method (with a constant gap 
parameter A) all these quasiparticle amplitudes are pos- 
itive, and hence all the resulting C matrices are positive 
definite, thus fulfilling the sufficient condition for the full 
diagonalization of C + C'=CC' . In fact, quasiparticle am- 
plitudes of different signs rarely occur in nuclear physics 
applications, cf. Ref. 0. This is so, because typical pair- 
ing forces couple the time-reversed states, and, in general, 
are always attractive. This shows that the Jordan struc- 
tures discussed here cannot be expected to be frequently 
encountered, and most often one will deal with the stan- 
dard canonical forms of Eqs. (^) and (|27]) , in which the 
only non-zero matrix elements are adjacent to the main 
diagonal. 

In summary, I have shown how to extend the results 
of Ref. in order to construct canonical basis in which 
two arbitrary complex antisymmetric matrices C and C 
acquire simultaneously canonical forms. This construc- 
tion completes the generalization of the classic Bloch- 
Messiah-Zumino theorem to the case of non-diagonal ma- 
trix elements calculated between fermion-pair conden- 
sates. 
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